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Abstract 

We introduce a fully frustrated XY model with nearest neighbor (nn) and 
next nearest neighbor (nnn) couplings which can be realized in Josephson 
junction arrays. We study the phase diagram for < x < 1 (x is the ratio 
between nnn and nn couplings). When x < 1/V2 an Ising and a Berezinskii- 
Kosterlitz-Thouless transitions are present. Both critical temperatures de- 
crease with increasing x. For X > I/V2 the array undergoes a sequence of 
two transitions. On raising the temperature first the two sublattices decouple 
from each other and then, at higher temperatures, each sublattice becomes 
disorderd. 
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A variety of two dimensional systems undergo a phase transition without a rigorous 
symmetry breaking, the Berezinskii-Kosterhtz-Thouless (BKT) transition [|l|. The transition 
is driven by the thermally excited vortices which form a two dimensional Coulomb gas 0. 
Josephson junction arrays are experimental realization of the two-dimensional XY model 
where the array's parameters can be modified in a controlled way. In the last decade there 
has been a great amount of work on the various aspects of the BKT transition in Josephson 
arrays 0] . Experimental studies are based on electrical resistance , two-coil inductance , 
and SQUID measurements. 

A magnetic field applied perpendicularly to the array leads to frustration If the 

flux piercing the elementary plaquette is half of the flux quantum $o = hc/2e, the system 
is called fully frustrated (FF) and undergoes two phase transitions related to the Z2 and 
f/(l) symmetries. The existence of the two critical temperatures T^'^ and T^/^^^ respectively 
has been extensively investigated both by analytical methods [HQ and Monte Carlo (MC) 
simulations |TD|-|T^. The complete scenario is not fully understood yet. There are numerical 
evidences either supporting the existence of two very close critical temperatures {T^'^ > 
T^^*^^^) with critical behavior typical of Ising and BKT transitions respectively |T3|-p!5| or the 



existence of a single transition with novel critical behavior [|TT|. 

Untill recently only the Josephson arrays with nearest neighbor (nn) couplings were 
studied. Recent theoretical and experimental |l7| works on infinite range array opened 



a new field of investigation in these systems. In this Letter we study the properties of a 
two dimensional FF Josephson array with both nn and next-to-nearest neighbor (nnn) cou- 
plings |18|. Proximity-junction arrays may be good candidates to experimentally probe the 



effects discussed in this work. They consist of superconducting islands in good electric con- 
tact with a metallic substrate. Due to the proximity effect there is a leakage of Cooper pairs 
into the normal substrate which extends over a temperature-dependent coherence length C,n 



{^N = fiVF/ksT or = ^JhD/2'KkBT for the ballistic and for the diffusive case respec- 
tively, vp is the Fermi velocity, D the diffusion constant and T the temperature). When 
^jv becomes comparable with the lattice constant of the array the nnn coupling becomes 
comparable with the nn coupling. Since is strongly dependent on the temperature, the 
nnn Josephson coupling may be observed cooling down the sample. The main results of this 
Letter are summarized in the phase diagram of Fig. ^. 
The system is defined by the Hamiltonian 

H = - E J^JCos{e,-eJ~A,^) (1) 

where Jij = J > for nn and Jij = xJ for nnn (x > 0), the symbol << ... >> refers to 
the sum over nn and nnn. For later convenience we introduce the gauge invariant phase 
difference (pij = 6i — 6j — A^. The variables 6i are the phases of the superconducting order 
parameter of the i-th island. The magnetic field enters through Aij = (27r/$o) // A-rfl (A 
is the vector potential). The relevant parameter which describes the magnetic frustration is 
f = J2 ^ijl (2vr), where the summation runs over the perimeter of the elementary plaquette. 
We study the case / = 1/2 on a square lattice. 

Ground States - The model of Eq. ( |l|) combines the characteristics of both the FF and 
unfrustrated XY models. While the elementary square plaquette is FF, the square plaquet- 
tes formed by nnn couplings are not frustrated. For x < Xq = 1/ \pl we find that the ground 
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state of this system is exactly the same as in the FF model without nnn couplings and 
is characterized by (ptj = ±7r/4 for all pairs of nn sites. This state combines continuous 
[?7(1)] and discrete (Z2) degeneracies. Its energy per site E = —^/2J is independent of nnn 
coupling xJ. Moreover if one considers a straight domain wall separating the two ground 
states with opposite orientations of chiralities it turns out that neither the form of such state 
nor its energy change with addition of nnn interactions. 

For X > xq the ground state is the same as in the absence of nn coupling when the 
system splits into two unfrustrated XF-models. The relative phase shift between the two 
sublattices can be arbitrary. The energy of this state E = —2xJ depends only on nnn 
coupling xJ and does not depend on nn coupling. 

At the special point x = xq the energies of both the above ground states coincide. 
Morover they can be transformed into each other by a continuous transformation without 
increasing the energy of the system. Therefore the manifold of the ground states also includes 
an additional set of eight-sublattice "intermediate" states which can be parametrised by a 
rotation angle X (x = ^ corresponding to low-x ground state and x = corresponding to 
high-x ground state with a particular relative phase shift between the sublattices) as it is 
shown in Fig. |l]. 

Phase Diagram - We studied the finite temperature behaviour of the model by means of 
the (low temperature) spin-wave free energy analysis and by Monte Carlo simulations. 

For X < Xo neither the spectrum of spin waves (in the long wavelength limit) nor the 
domain wall energy depend on x. Therefore we can expect only a weak dependence of T^^ 
and T,^(^) on x, due to the change of the effective interaction between the different types of 
fluctuations. 

For X > Xq the system (at finite temperatures) turnes out to be equivalent to two coupled 
XY models, the effective forth-order coupling between the two sublattices provided by the 
free energy of spin waves. Although this coupling is weak (always much smaller than the 
temperature) at low temperatures it is relevant and imposes the presence of a transition at 
T = Td. This transition separates the phases with coupled and decoupled sublattices. In 
the low temperature phase, where the two sublattices are locked, the spin-wave contribution 
to the free energy imposes a relative phase shift of ±7r/4 (or equivalently ±37r/4) between 
the two sublattices. At T > To a second phase transition of the BKT type takes place in 
each of the decoupled sublattices. For x 3> xq the temperature of this transition depends 
only on nnn coupling and it is proportional to x. 

The spin wave spectrum remains rigid down to x = Xq. This indicates (and it is confirmed 
by the MC simulations) that the critical temperature of this BKT transition remain finite 
when Below this temperature there is a transition between the FF low-x and 

the unfrustrated high-x phases. We evaluated numerically the spin wave free energy of the 
intermediated ground state as a function of x- This dependence is described by a convex 
function implying a first order phase trasition line separating the low-x and high-x quasi- 
ordered phases. 

The critical properties were investigated evaluating the helicity modulus F and the stag- 
gered chiral magnetization M by means of standard MC simulations for different x. The 
order parameter M, which controls the Ising-like transition, is defined as 

(2) 



rrii 
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where r^/a 
_ 1 

H 



+ sin I 



is the position vector 



in unit of lattice step a) of the site i and 
rrii = -^[sm (pij^ + sm (pjij2 + ^j2,j3 + 0i3,i) is the chirahty of the plaquette with center 
in {i^ + 1/2, iy + 1/2) and with site indexes i,ji,j2j3 (in clockwise order). 

The helicity modulus F = d'^T jdS^, used to signal the existence of a BKT transition, 
is defined through the increase of the free energy T due to a phase twist b imposed in one 
direction [^. In order to obtain a precise determination of T/^ and of the critical exponent 
V associated to the divergence of the correlation length we have calculated the Binder's 
cumulant of the staggered chiral magnetization M 



U 



M 



(3) 



Since U{T^'^,L) does not depend on lattice size L for large systems, T^^ can be identified 
without making any assumption on the critical exponents. Once a satisfactory estimation of 
T^'^ is obtained the critical exponent ly is estimated through a data collapsing with u left as 
the only free parameter. Estimation of Um have been obtained averaging, at least, 10^ • 
MC configurations by using a standard Metropolis algorithm. The largest lattice studied is 
L = 72. The result for x = 0.5 are shown in Figs.|,|a. We estimate fc^Tf ^ /J = 0.403±0.003. 
The data collapsing, shown in the inset of Fig.^ gives an estimate of l/u = 1.0 ± 0.1. The 
critical temperature T^^ decreases with increasing x for x < Xq. For x > Xq, there is no sign 
of a Ising-like transition (see Fig.^D). 

Following the procedure proposed in Refs. |]2T| , |22| , the critical temperature T^^(^) is esti- 
mated by using the following ansatz for the size dependence of F 



ttF 
27; 
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1 + 



1 



2(lnL-ln/o) 



(4) 



where Iq is a fit parameter. This critical scaling is based on the mapping between a neutral 
Coulomb gas and a XY model. Therefore Eq.(D can be used both as a test for the existence 
of a BKT transition and for a precise evaluation of the critical temperature. A very good 
scaling was obtained with 7 = 1 (the ordinary BKT transition) in the low-x phase and 7 = 2 
(corresponding to BKT transition on each of the two sublattices with lattice constant \/2) in 
the high-x phase. In Fig.^ we show this analysis for the 0.5 and x = l. Ty^^\x), 

as well as T^^{x), decreases with increasing x up to x ~ a:o. Our results cannot discriminate 
between the ^ = Tf and the Tf ^ > hypoth esis since the two temperatures are 

compatible within the numerical precision (the mean value of T/^ (x) remains always above 
the corresponding mean value of Tf (^)(x)). For x > 0.8, instead, Tj^^^^ (x) increases, quickly 
tending towards the value expected for x — *• cxd i.e. ksT^^^y {xJ) = 0.89. 

We finally discuss the transition related to the decoupling of sublattices in the high-x 
phase. The order parameter 5* 

can be chosen in the following gauge-invariant form: 



J2j,a ^j,j+ea can be defined on the bonds of the lattice and 



—ly'' exp 



i(-l)^'^exp \i{-iy=^+^y(f)^ 



This form of the order parameter is chosen in such way that for any of the high-x ground 
states the value of s will be the same for all the bonds. At low temperatures S manifests a 
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true long-range order. The MC simulations confirm the 7r/4 relative phase shift anticipated 
by the spin-wave analysis. By increasing the temperature the long-range order in S can 
be expected to disappear as a separate phase transition whereas the unbinding of vortex 
pairs in each of the sublattices has to occur at still higher temperatures. We performed MC 
simulations to evaluate the transition temperature T^. The results of this computation (not 
reported here) show that To is very close, but lower, than the transition to the disordered 
phase T^^^\x). More extensive simulations are needed to determine the critical behaviour 
of decoupling transition. The dotted line in Fig.^ shows the qualitative behaviour of the 
decoupling transition function of x. 

In conclusion we have introduced a frustrated XY model with nnn interaction. The 
model can be experimentally realized in Josephson junction arrays in a transverse magnetic 
field. Signatures of nnn Josephson couplings might already have been seen in specially 
designed setups [|3l. The analysis presented here leads to the phase diagram shown in 



Fig.^. For < X < Xq the critical temperatures T^^ and T^^^\x) decrease with increasing 
X. For X > Xq, there is no sign of a Ising-like transition and the system behaves like the 
unfrustrated XY model. At x ~ xq in the low T region there is a first order phase transition 
between the low and high-x phases. Finally for x > xq the array undergoes a sequence of 
two transitions. On raising the temperature, first the two sublattices become decoupled and 
then, at higher temperatures, each sublattice becomes disordered. 
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FIGURES 
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FIG. 1. The family of the additional (intermediate) ground states at the degeneracy point 
X = l/\/2 is characterized by the angle x- The value X = tt corresponds to low-x ground state 
while X = corresponds to high-x ground state with a particular relative phase shift between the 
two sublattices. The angle reported in the figure is the gauge invariant phase difference between 
two neighboring sites. 
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FIG. 2. The Binder's parameter Um vs. T for several lattice sizes L at x = 0.5. Errors are 
smaller than the symbols size. Excluding the data of smallest size (L = 24) one can estimate 
0.400 < ksT^'^/J < 0.406. Inset: Collapse of the data (excluding the L = 24 points). The scaling 
parameters are ksT^'^/J = 0.403 ± 0.003 and l/u = 1.0 ± 0.1. 
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FIG. 3. The staggered chiral magnetization M vs. T for a; = 0.5 a) and x = 0.8 b). For large 
L and low T, M goes to a non-zero value for x = 0.5 and vanishes for x = 0.8. The errors are 
smaller than the symbols size. The symbols for different lattices sizes are the same in a) and b). 
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FIG. 4. The size dependence of the hehcity modulus F for x = 0.5 and x = 1 at several 
temperatures T. The errors are smaller than the symbols size. The form of the plot [7T/ (TrT — 27r) 
vs. InL] is chosen in such a way that the scaling behaviour predicted by Eq. should correspond 
to a straight line with the slope equal to 1 (as shown by the over-imposed lines). The estimates 
for the critical temperature Tc^^^^ are 0.4 it 0.002 and 0.7 it 0.005 respectively. The error on the 
estimates is due to the used temperature mesh. 
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FIG. 5. The phase diagram for the XY model with nn+nnn interactions. Squares refer to T^'^, 
circles to Tc^^''; x is the ratio between nnn and nn couplings; J is the nn coupling. The errors are 
smaller than the symbols size (see text for details) . The dotted line shows the qualitative behaviour 
of the transition associated to the decoupling of the two nnn sublattices. 
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